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Abstract. We study the possibility to reproduce the statistical relations of the 
sunspot activity cycle, like the so-called Waldmeier relations, the cycle period - 
amplitude and the cycle rise rate - amplitude relations, by means of the mean field 
dynamo models with the fluctuating a-effect. The dynamo model includes the long- 
term fluctuations of the a-effect and two types of the nonlinear feedback of the mean- 
fleld on the a-effect including the algebraic quenching and the dynamic quenching 
due to the magnetic helicity generation. We found that the models are able to 
reproduce qualitatively and quantitatively the inclination and dispersion across the 
Waldmeier relations with the 20% fluctuations of the a-effect. The models with the 
dynamic quenching are in a better agreement with observations than the models with 
the algebraic a-quenching. We compare the statistical distributions of the modeled 
parameters, like the amplitude, period, the rise and decay rates of the sunspot cycles, 
with observations. 
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1. Introduction 

It is observed that the sunspot's activity is organized in time and latitude and forms 
the large scales patterns which are called the Maunder butterfly diagram. This pattern 
is believed to be produced by the large-scale toroidal magnetic field generated in the 
convection zone. Another component of the solar activity is represented by the global 
poloidal magnetic field extending outside the Sun and shaping the solar corona. Both 
components synchronously evolve as the solar 11-year cycle progresses. The global 
poloidal field reverses the sign in the polar regions near the time of maximum of sunspot 
activity. 

A remarkable feature of cyclic solar activity is that it is far to be just a cycle. 
Cycle amplitude and shape varies from one cycle to the other and prognostic abilities 
of any study of solar activity looks as its very attractive destination. Solar activity 
observations give various hints that various tracers of solar activity which are exploited 
to quantify the phenomenon demonstrate some relation one to the other what opens a 
possibility to predict future evolution of solar activity basing on available observations 
of other indices. Waldmeier jSS] pointed out at first this option (an inverse correlation 
between the length of the ascending phase of a cycle, or its "rise time", and the peak 
sunspot number of that cycle) and applied it, |37| . to give a prediction for the following 
cycle. The latter paper is in practice the first accessible (at least for German speaking 
readers) paper in the area. Later other relation of this this type was suggested and 
summarized as Waldmeier relations. This development was clearly summarized by 
and recently by |12] . The nature of the physical processes, that are manifested in 
the Waldmeier relations, is not clear, see discussion, e.g., in [HI HOI IE]- It seems to 
be remarkable, however, that these statistical properties of magnetic activity are also 
existed for the other tracers related with the sunspot activity (e.g., sunspot group and 
squares of sunspot groups, see IT^ |2]), and even for the other kind of the solar 
and stellar activity indices, e.g., for the Ca II index jSl]. The Waldmeier relations are 
considered as a valuable test of the dynamo models [T^ |S1 Ell- 

A natural way to push the understanding of the problem forward is to clarify the 
physics underlying Waldmeier relations. It is more or less accepted that cyclic solar 
activity is driven by a dynamo, i.e. a mechanism which transforms kinetic energy of 
hydrodynamical motions into magnetic one. Most of the current solar dynamo models 
suggest that the toroidal magnetic field that emerges on the surface and forms sunspots is 
generated near the bottom of the convection zone, in the tachocline or just beneath it in a 
convection overshoot layer (see, e.g., |32])- This kind of dynamo can be approximated by 
the Parker's surface dynamo waves |26]. The direction of the dynamo waves propagation 
is defined by the Parker- Yoshimura rule |38) . It states that for the aVL kind dynamo the 
waves propagates along iso-surfaces of the angular velocity. The propagation process 
can be modified by the turbulent transport (associated with the mean drift of magnetic 
activity in the turbulent media by means turbulent mechanisms), by the anisotropic 
turbulent diffusivity (see, [S]), and by meridional circulation [7J. A viewpoint, which 
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is an alternative to the Parker's surface dynamo waves is presented by the distributed 
dynamo with subsurface shear, e.g. [3j. The dynamo waves here propagates along the 
radius in the main part of the solar convection zone, [11]. The near surface activity is 
shaped by the subsurface shear. One more option is the flux-transport dynamo, e.g. 

[HE]. 

In the context of dynamo theory, the Waldmeier relations have to be explained by 
some mechanism which varies amplitude and shape of activity cycle and fluctuations a- 
effect are considered below as such mechanism. This idea extend the approach proposed 
in |29| to explain these relations by changing the magnitude of the a-effect. 

The physical idea underlying this mechanism can be presented as follows, a- 
coefficient is a mean quantity taken over ensemble of convective vortexes. Number 
N of the vortexes in solar convective shell is large however much smaller then, say, the 
Avogardo number, so fluctuations being proportional to N~^/'^ may be not negligible. 
Particular choice of N is obviously model dependent however if we take just for 
orientation = 10^ then A^~^/^ = 0.01. Taking into account that a is usually about 
1/10 of turbulent velocity we consider a dozen percent of a-fluctuations as a comfortable 
estimate. From the other hand, governing equations for large-scale solar magnetic field 
deal with spatial averaging and have to include a contribution of a-fluctuations, [13]. 

A straightforward application of the idea with vortex turnover time and vortex 
size as correlation time and length for a-fluctuations needs fluctuations much larger 
then mean a. [24j, |34j based on experiences in direct numerical simulations, e.g. [4J, 
and results of current helicity (related to a) observation in solar active regions, e.g. 
[39] considered a-fluctuations with correlation time comparable with cycle length and 
correlation length comparable with the extent of the latitudinal belts where sunspots 
occur to conclude that a reasonable a-noise of order of few dozen percents is sufficient 
to explain Grand minima of solar activity. The aim of this paper is to apply this idea 
to explain Waldmeier relations. 



2. Basic equations 



2.1. 2D model 

The dynamo model is based on the standard mean-field induction equation in perfectly 
conductive media [19j: 

= V X (5+U X B) 

at 



where £ = u x b is the mean electromotive force, with u, b being the turbulent 
fluctuating velocity and magnetic field respectively; U is the mean velocity (differential 
rotation). The axisymmetric magnetic filed: 

r smt7 
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Figure 1. Parameters of the solar convection zone: a) the contours of the constant 
angular velocity plotted for the levels (0.75 — 1.05)rio with a step of 0.025ilo; i^o = 
2.86- 10~''s~^; b) turnover convection time Tc, and the RMS convective velocity u'^ and 
the background turbulent diffusivity -q"^^ profiles; c) the radial profiles of the a-effect 
tensor components. 



9 - polar angle. We have used the expression for S obtained by (hereafter P08) and 
write it as follows: 

£i = (aij + 'jij) Bj - TjijkVjBk. (1) 

Tensor represents the alpha effect, including the hydrodynamic and magnetic helicity 
contributions, 

ci^, = (1 + i) sin^ ^ajf ) + af\ (2) 

where the hydrodynamical part of the a-effect, ai^ , C is the noise, and the quenching 
function, if) a-, are given in Appendix (see also in |28]). The hydrodynamic a-effect term 
is multiplied by sin^ 6 {6 is co-latitude) to prevent the turbulent generation of magnetic 
field at the poles. The contribution of the small-scale magnetic helicity x = ^^'b (a is a 
fluctuating vector-potential of magnetic field) to the a-effect is defined as a]^^^ = C-^^Xi 
where coefficient C^^^ depends on the turbulent properties and rotation, and is given 
in Appendix. The other parts of Eq.([T]) represent the effects of turbulent pumping, 7jj, 
and turbulent diffusion, rjijk. They are the same as in PKll. We describe them in 
Appendix. 

The nonlinear feedback of the large-scale magnetic field to the a-effect is described 
as a combination of an "algebraic" quenching by function ipa (P) (see Appendix and 
\29{). and a dynamical quenching due to the magnetic helicity conservation constraint. 
The magnetic helicity, x 5 subject to a conservation law, is described by the following 
equation [181 El 133] : 

where Tc is a typical convection turnover time. Parameter controls the helicity 
dissipation rate without specifying the nature of the loss. It seems to be reasonable that 
the helicity dissipation is most efficient in the near surface layers because of the strong 
decrease of Tc (see Figure lb). The last term in Eq.([3]) describes the diffusive flux of 
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Figure 2. The typical time-latitude and the time-radius (at the 30° latitude) diagrams 
of the toroidal field (grey scale), the radial field (contours at left panel) and the poloidal 
magnetic field (contours at the right panel) evolution in 2D1 model (see Table 1). The 
toroidal field averaged over over the subsurface layers in the range of 0.9 — O.99i?0 , 
the radial field is taken at the top of the convection zone. 

magnetic helicity j20]. We use the solar convection zone model computed by |32| . in 
which the mixing-length is defined as £ = aMLT |A*^p)| , where A^^-* = Vlogp is the 
pressure variation scale, and Umlt = 2. The turbulent diffusivity is parametrized 

in the form, rjT = Cr^rl^\ where rf^'' = — is the characteristic mixing-length 

3 

turbulent diffusivity, £ and u' are the typical correlation length and RMS convective 
velocity of turbulent flows, respectively and C,, is a constant to control the intensity 
of turbulent mixing. In the paper we use = 0.05. The differential rotation 
proflle, Q = Qofn{x, fi), x = t/Rq, fi = cos6 is a modifled version of an analytical 
approximation to helioseismology data, proposed by j2], see Fig. la. 

We use the standard boundary conditions to match the potential field outside and 
the perfect conductivity at the bottom boundary. As discussed above, the penetration 
of the toroidal magnetic field in to the near surface layers is controlled by the turbulent 
diffusivity and pumping effect. For magnetic helicity, similar to [H] and [2T], we use 
the time dependent conditions provided be EqjS] and the helicity flux conservation the 
condition = is applied at the bottom and at the top of domain. The latter gives 
a smooth transfer for solutions with and without the diffusive helicity flux. 

The left panel on the Fig. [2] shows the typical the time-latitude diagram for the 
toroidal magnetic averaged over the subsurface layers 0.9 — 0.99/?© and the radial 
magnetic at the top of the integration domain. The right panel shows the time-radius 
the time-radius diagram for the toroidal an poloidal magnetic fleld evolution at 30° 
latitude. 

We demonstrate it by Fig. 3 which shows the time-latitude diagrams for toroidal 
and radial magnetic fleld evolution for the models IDl, 1D3 and 2D1. For the latter 
model we show the toroidal magnetic averaged over the subsurface layers 0.9 — O.99i?0 
and the radial magnetic fleld is given for the top of the integration domain. For the 
model 2D1 we show the time-radius diagram for the toroidal an poloidal magnetic fleld 
evolution at 30° latitude. The other models listed in Table 1, having the same general 
patterns of the magnetic fleld evolution, are differed from the models shown on the 
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Fig. 3 in some details (mostly associated with magnetic helicity evolution). 

2.2. ID model 

For comparison with the previous studies and also to study how the additional dimension 
affect the statistical properties of the dynamo we consider the ID model similar to that 
studied by [21]: 

BA (9/1 dA\ 

— = sin ^ ((1 + cos ei,^{B) + X) 5 + sin ^- (^ — — j - r^czA, (4) 

where the large-scale radial shear VL{9) = dfl/dr. The ID model employs two 
possibilities for the shear profile. In one case we put Cl{6) = 1, that give us the model 
explored by ^3]. In another case we use 

n=^{5sm'e-A), (6) 

which is suggested by [16]. In agreement with the helioseismology results for the bottom 
of the convection zone, this profile is positive in equatorial regions and negative near the 
poles. The magnetic field strength in Eq.(|5]) is measured in the units of the equipartition 
magnetic field strength and the time is normalized to the typical diffusive time, R^/ifjp . 
The evolution of the magnetic helicity for the ID model is governed by equation: 

2 dA d sin 6B x Vx ^ f ^ 

^ sin^ e ~de d~e ^ ^ '^de 

In what follows we will discuss the ID models with the constant shear, because they 
are more relevant to compare with observations. The differences in results for the ID 
models with the variable shear given by Eq.(|6]) will be briefly mentioned in subsequent 
sections. 

Summarizing, we exploit much more detailed and realistic dynamo models then 
|24j . |34] . Our point is that Waldmeier relations are a much more delicate phenomena 
rather Grand minima and the bulk of our knowledge concerning recent solar cycles is 
much more rich then that one for remote past when Grand minima took place. 

2.3. Noise model 



The noise, ^, contributes in the hydrodynamic part of the a-effect (see, Eqs.( 2|[4 )) 



Following to [23] the models employ the long-term Gaussian fluctuating ^ of the small 
amplitude with RMS deviation given in the Table 1 (last column). The time of the 
renewal of the ^ is equal to the period of the model. The random numbers were 
generated with help of the standard F90 subroutine quality of contemporary standard 
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Figure 3. The left panel shows the time-latitude diagrams of the toroidal field (grey 
scale) and the radial field (contours) for the IDl model (see Table 1). The right panel 
shows the estimated sunspot number in the the separated cycles in the IDl model (see 



Section 2.4) 
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Cw 
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IDl 


1 








10 


3 


1200 


0.15 


2D1 




Eq.(2) 


10-5 


200 


800 


1 


0.15 


2D2 




Eq.(2) 


0.3 


10*^ 


200 


1 


0.15 



Table 1. Parameters the dynamo models: the type of the nonlinear quenching of the 
a-effect, if the magnetic helicity is x = then the model employ only the algebraic 
quenching which is described by ^a&Tid otherwise by the dynamic quenching due 
to magnetic helicity described by Eq.([3]) or Eq.([7]); rj^/rjT is the ratio between the 
turbulent magnetic helicity diffusivity and the turbulent magnetic diffusivities; the 
profile of the shear in the ID models; the a-effect parameter in the 2D models; the 
parameter i?,^. controls the helicity dissipation rate; the parameter Bo/Beq controls 
the sunspot number parameter in the ID models. It is the ratio between the typical 
strength of the toroidal magnetic field producing the sunspots and the equipartition 
magnetic field strength; Bq is the typical strength of the toroidal magnetic field 
controlling the sunspots number parameter in the 2D models; Cw is the parameter 
to calibrate the modeled sunspot number relative to observations; a is the standard 
deviation of the Gaussian noise in the model 



noise generator subroutine is shown to be sufficient for such kind of modelhng, see e.g. 
[1]). It would be more realistic to consider the renewal time as the fluctuating quantity 
as well, but we would like to separate this effect for the different study. Also, we found 
that the models which employ the magnetic helicity effect show the very intermittent 
long term behaviour. This makes the analysis procedure (e.g., division to subsequent 
cycles) more complicated. We isolate ourselves from these phenomena by considering 
the noise models with the lower RMS in case if the magnetic helicity is employed. [|] 

I In part, the given problem is likely due to the very rough model for the Wolf number, see Eq.(p|). 
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2.4- The sunspot cycle model and the Waldmeier relations 

In the paper we define the Waldmeier relations as the set of the mean properties of the 
sunspot cycle. We will deal with the following properties of the Wolf sunspot number 
(which is taken either from observational database or simulated from the model): the 
relation between period and amplitude of the same cycle, the relation between rise rate 
and amplitude of the cycle and the shape of the sunspot cycle, characterized by the ratio 
between the decay rate and the rise rate in the cycle. The other kind of relations, like the 
link between the rise time and amplitude of the cycle, can be considered as the derivative 
from the above relation. For comparison with other analysis of the observational data 
and also with the results of the dynamo models presented by Karak and Choudhuri[8] 
we show the results for the rise time of the cycles as well, the relation between the 
rise time and amplitude of the cycle and the relation between the cycle amplitude and 
period of the preceding cycle (see, j35] and [I2|)- The amplitude of the cycles is defined 
by difference between the maximum sunspot number and the sunspot number in the 
preceding minimum. Even for the harmonic cycles the latter differs from zero due to 
the spatial overlap in subsequent cycles. The period of the cycle is equal to the time 
between the subsequent minima, the rise time of the cycle is defined by the difference 
between the moment of the cycle maximum and the moment of the preceding minimum 
of the cycle. The rise rate is defined as the ratio between the difference of the sunspot 
number amplitude during maximum and minimum of the cycle and the rise time of the 
cycle. The similar definition is for the decay rate of the cycle. 

Remind that sunspots are not directly presented in dynamo models and we have to 
relate its number to a quantity involved in a dynamo model under consideration. We 
assume that the sunspots are produced from the toroidal magnetic fields by means of 
the nonlinear instability and avoid to consider the instability in details. To model the 
sunspot number W produced by the dynamo we use the following anzatz 



where for the 2D models (-Bmax)5i is the maximum of the toroidal magnetic field strength 
over latitudes averaged over the subsurface layers in the range of 0.9 — O.99i?0 and for 
the ID models (i?max)5L is simply the maximum of the toroidal magnetic field strength 
over latitudes; Bq is the typical strength of the toroidal magnetic field that is enough 
to produce the sunspot; Cw is the parameter to calibrate the modeled sunspot number 
relative to observations. The all parameters which were employed in the different models 
are listed in the Table 1. 

In the dynamo models we explore the effect of the Gaussian fluctuations of the 
a-effect, or parameter Cq with the typical time equal to the period of the cycle and the 
standard deviations less than 0.2Cq. In the models presented here we fix the standard 
deviation to O.ISCq. 

For comparison with simulation we use the smoothed data set from |30| which starts 
at 1750. Choosing this data set we appreciate that in principle Waldmeier relations can 
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-82.1X+617.4 


-25.6X+207.5 


-33.0X+252.8 


-26.7X+234. 


-16.1X+165.4 


Amplitude 


±18.3 


±35.3 


±22.7 


±24. 


+38.5 




-0.84 


-0.49 


-0.47 


-0.75 


-0.48 


Rise Rate - 


1.0x+4.0±3.1 


0.43x+3.3±2.2 


0.68X+1.6+1.6 
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0.9 


0.92 
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Table 2. First five rows contain information for the mean and variance (standard 
deviation) for the parameters of the sunspot cycles in the different data set. The shape 
of the cycle is defined as ratio between the decay rate and the rise rate of the cycle. 
Last five rows show the linear fits with the mean-square error bar and the correlation 
coefficient. In the relation Period- Amplitude (a) we compare the cycle amplitudes to 
period of the preceding cycle (see |12l I35| l. and in the relation Period-Amplitude (b) 
we compare these parameters for the same cycle. 

be valid for normal cycle only and their applicability to epochs of Grand minima of solar 
activity must be addressed separately. Available instrumental data concerning solar 
activity in XVII - early XVIII centuries gives a limited possibility only to address this 
important point which obviously is out of the scope of this paper. From the other hand, 
there are various indirect (mainly isotopic) tracers of solar activity which give a limited 
information concerning its shape over much longer time interval rather instrumental 
data. Our point is that Waldmeier relations and the regularities of such long-term 
time series (see, e.g..|23[ [22] ) have to be discussed in a separate paper and here use 
as an illustrative example the extended time series of the sunspot data proposed by 
|25| (referred hereafter as NIMV). These data sets are shown on Fig. |4| The Table 2 
contains the linear fits and correlations between the different parameters of the cycles for 
observational data sets and for the dynamo models as well. In particular, the parameters 
of the relation between rise time and amplitude and parameters of the Amplitude-Period 
effect (a) and (b) (associated with period of the preceding and the same cycle) for SIDC 
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Figure 4. The sunspot data sets. Upper raw: left - SIDC and right - [1S](NIMV), 
lower raw - corresponding cycles distributions 

data set are in a good agreement with the results by Vitinskij et al |35| and Hathaway 
et al[T2]. The similar conclusion can be done if we compare our analysis for SIDC data 
set for the relation between rise rate and amplitude of the cycles with analysis given by 
Vitinskij et al [35j . 

3. Results 

The typical time-latitude diagrams for the dynamo models were shown in Figures |2] and 
|3j The shape of the simulated sunspot cycles in IDl model can be seen on the right 
panel Figure |3j The simulated sunspot cycles for the 2D1 and 2D2 models are shown on 
the the Figure [s] We can conclude that the shape of the simulated sunspot cycles (and, 
perhaps, the associated Waldmeier relations) is directly related with the spatial shape 
of the toroidal magnetic field evaluational patterns. For example, in the IDl model 
the maximum of the butterfly diagram is very close to equator and butterfly wing is 
elongated toward the pole. In such a pattern of the toroidal magnetic field evolution the 
decay phase of the sunspot activity is shorter than the rise phase. The opposite situation 
is in the models 2D1 and 2D2. The physical mechanisms which produce the short rise 
and the long decay of the toroidal magnetic field activity were discussed recently by 
Pipin and Kosovichev [29j. 

To proceed further we would like to discuss the statistical properties of the cycle 
parameters those involved in the Waldmeier relations. The ID models have the much 
less cycle period than diffusive time of the system. Therefore, we scale the periods of 
these models by factor ~ 50 . The Table 2 show the results for the mean and the 
variance (standard deviations) for the period, amplitude, rise rate and the shape of the 



The fluctuating a- effect and Waldmeier relations in the nonlinear dynamo models 11 




'^O 2 4 5 8 10 12 14 16 
YR 







'0 2 4 6 8 10 12 14 
YR 



16 



Figure 5. Left panel shows cycles distributions for the model 2D1 and the right panel 
- model 2D 2. 



sunspot cycles in the different data sets. From that Table we see that the IDl model 
has the smaller variance in the period, amplitude and rise rate of the cycles as compared 
to the others data sets. The shape asymmetry of the cycles in IDl is opposite to the 
others cases as well. Also we can see that the mechanism of the helicity loss in the 
dynamo model influences the mean and variance of the sunspot cycles parameters. In 
particular, the model 2D2 with the increased diffusive loss of the magnetic helicity has 
the lower variance of the period and amplitude of the sunspot cycles and has the more 
symmetric shape of the cycle as compared to the model 2D1. The difference in the 
synthetic data set of the sunspot cycles provided by NIMV as compared with the SIDC 
is likely due to the fact that the SIDC data set does not cover the periods with low 
magnetic activity. This argument is also applied if we compared NIMV and, e.g., 2D1 
model. The parameters of the 2D1 model does not allow to have the extended periods 
of time with very low sunspot cycles. 

The difference of the the statistical properties of the given data set can be seen in 
further detail using the cumulative distribution probability functions. The cumulative 
distributions are constructed as follows. At the beginning, we sort each distribution for 
each parameter and each model in increasing order. After this we compute the following 



where Pi is the parameter under consideration (say, the cycle period) having the order 
number i (after sorting the set in increasing order) and is the total number of the 
instances of the given parameter in the set. Equation ^ approximate the probability 
for the parameter P to have the values in interval between Pmin and Pj. The accuracy 
of the approximation improves under N oo. We will use the log-normal cumulative 
distribution constructed on the base of the SIDC data set as the reference distribution. 
The SIDC data set has only 23 instances of the sunspot cycles. To construct the reference 
log-normal distribution we use the standard mean and variance of the cycles parameters 



CDF(P,) 



(9) 



N 
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Figure 6. CDF distributions, red line - SIDC the data set, blue line - the data set 
from [25j. 

(period, amplitude, rise rate and asymmetry) given in the Table 2. Then take the 
natural logarithm of them and construct the log-normal distribution of the length 1000 
using those mean and variance. The results are shown on the Fig. 6. 

It is clearly seen that log-normal distribution is a good fit for the distributions of the 
sunspot cycles period in the SIDC data set and also for model 2D2. The difference of the 
SIDC data set from the log-normal distribution is seen in the probabilities distributions 
for the rise rates and the shape of the cycles. It is, however, unclear if these differences 
is due to the limited data set of cycles covered by SIDC. The data set produced by 
the models and the NIMV data set can be equally well approximated by the log-normal 
distributions (with different mean and variance). For the dynamo models, the difference 
between the distributions computed by Eq.(|9]) and the log-normal approximations for 
them is less visible than for SIDC and NIMV sets. 

Fig. [7] shows the Waldmeier relations for the IDl and 2D1 models together with 
their linear fits and also fits for the SIDC and NIMV data sets. The parameters of 
the linear fits are summarised in the Table 2. It is seen that the model 2D1 is well to 
reproduce the SIDC data set, and the difference to the NIMV data is not very large. 
The correspondence of the 2D2 model to the SIDC and the NIMV is not as good as 
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Figure 7. The Waldmeier relations for IDl (left) and 2D1 (right) models. The linear 
fits are shown the solid lines, the dashed lines shows the fits for the SIDC data and 
the dash-dot hne - for the NIMV data. 



for the 2D1 model. This is also can be expected by results presented in Fig. [6] and by 
Table 2. Finally, we can conclude that IDl model has only qualitative agreement for the 
relations between the rise rate - amplitude, and the period - amplitude of the sunspot 
cycles. 

4. Discussion and Conclusions 

In the paper we have studied the possibility to reproduce the statistical relations of the 
sunspot activity cycle, like the so-called Waldmeier relations, by means of the mean field 
dynamo model with the fluctuating a-effect. The dynamo model includes the long-term 
fluctuations of the a-effect. The dynamo models employ two types of the nonlinear 
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feedback of the mean-field on the a-effect including the algebraic quenching and the 
dynamic quenching due to the magnetic helicity generation. The paper presents the 
results for three particular dynamo models. 

The presented ID model is similar to model discussed by Moss et al. p4j. It uses 
the constant shear and the algebraic quenching of the a-effect. The results for this model 
disagree with observations (SIDC data set) about the shape of the simulated sunspot 
(decay rate is higher than rise rate) even though it is qualitatively reproduce the basic 
Waldmeier relations for the Rise Rate-amplitude and the cycle Period-amplitude (see 
left column in Fig. [?]). It was found that the variance of the cycle parameters in the 
long-term evolution is less than in 2D models. It is interesting, that under the level 
of noise the ID models involving the magnetic helicity show the smaller mean even 
though having the stronger variances of the simulated sunspot parameters. Although 
we could scale the mean parameters of those models to the observational values, we 
did not present the results for these models because they have the Waldmeier relations 
which are quantitatively the same as those presented for IDl model in Table 2 and 
Fig. [71 

We checked the ID models with the spatially variable shear like that suggested by 
Kitchatinov et al. |16J. In agreement with the helioseismology results, the given ID 
models have the realistic latitudinal profile of the shear (see Eq.(|6])). Although, these 
models qualitatively reproduce the relation between the rise rate and amplitude of the 
cycle, they fail with the other kind of relations, having the positive correlation between 
the period and amplitude of the cycle and the equal rate for the rise and decay phase 
of the simulated sunspot cycles. 

Similar to the ID cases the magnetic helicity contribution to the a-effect results to 
decrease of the toroidal magnetic field strength and to growth the variance of the cycle 
parameters in the long-term evolution of the magnetic activity. The strong variance of 
the cycle parameters is expected from SIDC data set and from NIMV as well. For this 
reason in the paper we discuss the 2D model which involves the effect of the magnetic 
helicity. The 2D models employ two different description for the magnetic helicity loss, 
to overcome the problem of the a-effect catastrophic quenching. The term — x/i?^rc in 
Eq.([3| describes the magnetic helicity loss with the dissipation rate {tcR^)~^ without 
specifying the nature of the loss. Note, that Tc is varied from about 2 months near bottom 
of the convection zone to a few hours at the top of the integration domain (which is 
O.99i?0). Thus, for the = 200 used in the model 2D1, the typical decay time for the 
magnetic helicity is varied from about 4 solar cycles at the bottom of the convection 
zone to a time which is less than one month at the top of the convection zone. It is not 
clear if this simple description is satisfactory approximation for the magnetic helicity 
loss. Therefore we checked the alternative possibility using the diffusive helicity flux. 
Although, the model that employ the diffusive helicity flux is in satisfactory agreement 
with SIDC data, the correspondence to observation in this model is not as good as for 
the model 2D1. We find the the variance of the cycle parameters in the model 2D2 is 
less than in the model 2D1 while the SIDC and NIMV data sets show higher variances 
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than the model 2D1. 

The detailed comparison the results of our models with those given by Karak 
and Choudhuri |8] is not possible, because we have used a different definition for the 
amplitude of the cycle and the rise time. They did not give the results for the linear 
fits coefficients and only provide the correlation coefficients in the Waldmeier relations 
involving the Rise Rate-amplitude and the Rise Time-Amplitude of the cycle. Bearing 
in mind the differences in definition that their "high diffusivity model" with fluctuating 
meridional circulation is comparable with our 2D1 and 2D2 models. It is not clear 
however what is the typical shape of the cycle in their model. This is an important 
issue as we have seen in example given by model IDl. It has qualitative agreement 
with SIDC data about the period - amplitude and the rise rate - amplitude relations 
even-though having the rise time of the cycle greater than the decay time. 

In the models under consideration, the asymmetry between the ascent and decent 
phase of the sunspot cycle is inherent from the pattern of the toroidal magnetic field 
activity. In particular, the IDl model has the toroidal magnetic field butterfly diagram 
with maximum located very close to equator. Therefore, applying the definition Eq.(|8]) 
for this type of the toroidal magnetic field evolutional pattern we obtain the decent 
phase of the sunspot activity shorter than the ascent phase. The opposite situation 
is in 2D models. There, we relate the sunspot activity with the toroidal field in the 
subsurface layers. The turbulent diffusivity in the model decrease outward this leads 
to increases the decay time when the toroidal field gets closer to the surface (see |29]). 
We find that the effect of the magnetic helicity on the a-effect can amplify or saturate 
the asymmetry of the cycle shape depending on the mechanisms of the helicity loss 
employed in the model. 

It is expected that the nonlinear dynamo mechanisms affect both the magnetic cycle 
profile and the statistical properties of the cycles. The paper illustrates the impact of 
the non-linear a-effect for the algebraic and the dynamic non-linearities. Recently, 
Kitchatinov & Olemskoy |17| suggested that the non-linear diffusion could promotes the 
events similar to the Maunder minimum provided there are the small fluctuations in the 
a-effect. This mechanism does not work in our models, because on the rise phase of the 
cycle, the growing toroidal magnetic fields results to the turbulent diffusivity quenching 
and this effect makes the rise phase of the cycle longer, i.e., the smaller turbulent 
diffusion, the longer evolutionary time scale. The opposite situation is expected for the 
decay phase of the magnetic cycle. 

The comparison of the SIDC data set and the synthetic data set provided by 
Nagovizyn et al. |25|(NIMV) reveals the significant difference in the statistical properties 
of the cycle parameters. This seems to be a result of the wider cycle variations range 
covering by the NIMV data set. The model presented in the paper don't cover the 
variations seen in NIMV because the selected models almost have no the extended 
events with low cycles like the so-called Maunder minimum which were observed during 
the 16-th century. This motivated us to extend our study and explore the models which 
have the more intermittent variations of the sunspot cycle. This work is planned for the 
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future papers. 

Summarizing the main findings of the paper we conclude as follows. We found 
that the dynamo models, having the reasonably good the time-latitude diagram of the 
toroidal magnetic field evolution, are able to reproduce qualitatively the inclination 
and dispersion across the Waldmeier relations with less than 20% Gaussian fluctuations 
of the a-effect. The 2D models have better agreement with observations than ID. In 
particular, ID models fail to reproduce the asymmetric shape of the sunspot cycle with 
short rise and long decay phases. The statistical distributions of the cycle parameters 
show the log-normal probability distributions for the all data sets analysed in the paper. 
The parameters of these distributions are different for all data sets. Again the ID 
model is significantly different from others in this sense. The 2D model that employs 
the simplest form of the helicity loss via the term — x/i^^r^ agrees well with the SIDC, 
even-though the long-term variations in this model is not intermittent enough, and 
this seems to be a reason for its difference to the NIMV data set in some aspects. 
The employ of the diffusive loss in the magnetic hehcity evolution equation results to 
decreasing in the variations of the cycle parameters. The further study of the magnetic 
helicity transport mechanisms should clarify the likely candidates which are responsible 
for the magnetic hehcity loss from the dynamo region. We have seen that the analysis 
of the statistical relations of the sunspot cycle may provide the valuable diagnostic tool 
for this study. 
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5. Appendix 

We describe some parts of the mean-electromotive force. The basic formulation is given 
in P08. For this paper we reformulate tensor a - ^ , which represents the hydrodjTiamical 
part of the a-effect, by using Eq.(23) from P08 in the following form, 



a - ■ 



= 5., {3VT (/ff (e • A(-)) + fif (e ■ A^")) ) } + (10) 
+ e.e, {3,T (/^^ (e-A(^))+/r (e-A^")))} 



The contribution of magnetic helicity x = ^i-b (a is a fluctuating vector magnetic field 
potential) to the a-effect is defined as = C^ifxi where 

Cff = 2/i'^)5,,^ - 2/(")e.e,^. (11) 



The turbulent pumping, 7j is also part of the mean electromotive force in Eq.(23)(P08). 
Here we rewrite it in a more traditional form (cf, e.g., ), 

= ^VT {ft^Al'^ + fi"^ (e ■ A(^)) e„} - S^Tfi'^e^e.^^e^M^l (12) 
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The effect of turbulent diffusivity, which is anisotropic due to the Coriohs force, is given 
by: 

mjk = 3?7t I (^^ft^ - fi'^^^ Sijk - '^/["^eiCnenjk'^ ■ (13) 

Functions f^'^'^l^ depend on the Coriolis number fi* = 2tc^o and the typical convective 
turnover time in the mixing-length approximation: = i/u'. They can be found in P08. 

The turbulent diffusivity is parametrized in the form, rjjn = Cr^ril^\ where r^^^ = — is 

3 

the characteristic mixing-length turbulent diffusivity, u' is the RMS convective velocity, 
I is the mixing length, is a constant to control the intensity of turbulent mixing. 
The others quantities in Eqs. ( |l0||l2||l3| are: A{p) = Vlog p is the density stratification 



scale, A(") = Vlogf ?7y''') is the scale of turbulent diffusivity, e = r2/|fi| is a unit 
vector along the axis of rotation. Equations ( [I0|l2|13[ ) take into account the influence 
of the fluctuating small-scale magnetic flelds, which can be present in the background 
turbulence and stem from the small-scale dynamo (see discussions in). In our paper, 

the parameter e = — which measures the ratio between the magnetic and kinetic 
/iopu2 

energies of fluctuations in the background turbulence, is assumed equal to 1. This 
corresponds to the energy equipartition. The quenching function of the hydrodynamical 
part of a-effect is deflned by 

Note, in notation of P08 ipa = ~3/40g°'*, and /3 



Lb 



